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ABSTRACT 



The Belle experiment at the KEK B factory has collected 93 fb^^ 
of electron positron collisions at ^/s = 10.6 GeV. This has produced a 
sample of 85 million BB meson pairs that can be used to study CP 
violation in rare (and not so rare) B decay modes. Here I report on 
a measurement of indirect CP violation in the decay B^ J/ipKg, 
as well as time dependent CP asymmetries in rarer modes such as 

— > TT^TT^ jTj'Kg and 4>Kg. I summarise the prospects for improving 
the precision on these and related measurements. 
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1 CP Violation in B Decay 



When CP violation was first observed in neutral kaon decay, in the early 1960s, 
it shook the foundations of particle physics. It had previously been assumed that 
the combination of charge conjugation and a parity transformation left all known 
particle interactions invariant, despite the fact that weak interactions violated 
parity alone. Over the following three decades CP violation in the system 
was measured with ever increasing precision in an attempt to pin down its source. 
In the 1970s Kobyashi and Maskawa0 showed that the Standard Model could 
accommodate CP violation in a three quark weak mixing matrix, Vckm- The 
single non-trivial phase in such a 3 x 3 unitary matrix could explain the small 
effect first seen in meson decay, where CP violation was observed at the 10"'^ 
level. It was suggested,! in the early 1980s, that the comparable amplitudes for 
the direct decay of mesons into CP eigenstates and the mixing of B^B^ mesons 
would make neutral B meson decay an ideal place to observe large indirect CP 
violating effects. 

One way to understand the magnitude of CP violation predicted by the CKM 
model, in neutral B meson decay, is to consider the unitarity relation between the 
first and third columns of Vckm: 

v:,vu + v:,v,, + v:,Vuci = o. 

While each term in this expression is relatively small ((9(sin^ ^c)), they are all the 
same size. When plotted in the complex plane (see fig. ^ one expects significant 
angles at each apex because the sides of the triangle have similar lengths. In neu- 
tral kaon decay the corresponding unitarity triangle has two sides that are much 
larger than the third - making the decay rates for mesons much larger, but 
making the angles, and hence the observable phases small and more challenging 
to measure. 

In B meson decay CP violating phases are most readily observed through the 
indirect mixing of two amplitudes. Given a CP eigenstate accessible to both B^ 
and -B" decays - such as J/ipKg - one can observe the interference between the 
direct decay amplitude for: 

B'{B') J/^Kl- 
and the amplitude for the same decay preceded by B^ mixing: 
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Figure 1: Unitarity of the CKM matrix relates its different elements. The study 
of CP violation in B meson decay involves the CKM elements depicted in this 
triangular relationship. 



When one includes the relative phase between these two amplitudes we get an 
expression for the time dependent CP asymmetry: 



AcpiAt) = 



- C,K sin 201 sin AmAt. [1] 



Where the asymmetry in the decay rate between and B^ mesons is proportional 
to the B^ mixing rate (sin AmAt), the CP eigenvalue of the final state, C,k, 
{^x = ~1 for J/ipKg and +1 for J/ipK^) and sin20i, the angle at the lower right 
apex of the unitarity triangle shown in fig. |T]. 



2 The KEK-B Collider and Belle Detector 

The main experimental challenge in measuring CP violation in B^ meson decay 
lies in the fact that B^ mesons decay much more quickly than mesons, hav- 
ing proper flight distances of fractions of a millimeter. Furthermore the decay 
rate to experimentally accessible CP eigenstates are much smaller, on the or- 
der of (9(10~^) for B^ mesons; while essentially 100% of meson decays are 
to identifiable CP eigenstates. While the very much larger CP violation in B^ 
decay makes up for some of this it remained a significant experimental challenge 
to observe CP violation in B^ meson decay. 

The key ingredient is having a high luminosity source of B mesons - a S 
factory. Two dedicated machines were built in the late 1990s to address this. The 




Figure 2: The Belle experiment has been collecting data over the last two years. 
This plot shows the integral luminosity over that period (lower line) compared 
to the luminosity collected by the BaBar experiment at PEPII (upper line). The 
data samples are similar at this point. 



KEK-B e^e~ accelerator complex collides beams of 8 GeV electrons and 3.5 GeV 
positrons, to produce T(4S') mesons at y/s = 10.6 GeV. These, in turn, decay into 
BB meson pairs. Using an 11 mrad crossing angle KEK-B is able to collide beams 
with currents in excess of 1 Ampere, with tolerable backgrounds and luminosities 
approaching 10'^^ per cm^ each second. An integrated luminosity comparison 
between the KEK-B factory and its competitor, PEP-II at SLAG, is shown in 
fig. ^ While the two machines have delivered similar integrated luminosities since 
they came online two years ago, the recent instantaneous luminosities (the slope 
of the curve) in the KEK machine bodes well for upcoming data-taking. 

The colliding beam energies are asymmetric so that the T{4S) meson produced 
is boosted along the beam direction. This results in the daughter B mesons also 
being boosted, with /?7 ^ 0.45, separating the two B meson decay vertices in the 
detector. The boost introduces some asymmetry in the design of the experiment. 
A cross-section of the cylindrical Belle detector is shown in fig. ^. The detector 
elements are arranged asymmetrically around the interaction point to increase the 
acceptance for B meson decay products. 

Working outwards from the collision point B meson decay products encounter 
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Figure 3: An elevation view of the Belle experiment. 



a three layer silicon vertex detector (SVD) that measures charged particle trajecto- 
ries with 55 /im precision (at 1 GeV/c) and the separation of h decay vertices with 
a precision of 100 /xm in z. They next pass through a Helium filled drift chamber 
(CDC) that measures track momenta with a precision of ap/p = {0.2p © 0.3)% (p 
measured in GeV/c) providing excellent mass resolution for B decays into charged 
daughter tracks. This is followed by a Cesium-Iodide crystal calorimeter (Csl) that 
has better than 2% energy resolution for 1 GeV photons. Belle's particle identifi- 
cation system includes an aero-gel Cerenkov counter system (PID) and a time of 
flight system (TOF) that can distinguish kaons from pions up to 3.5 GeV/c with 
90% efficiency and fake rates of less than 5%. These systems are followed by the 
solenoid coil and then a and muon detection system (KLM) that identifies 
muons with less than 2% fake rate above 1 GeV/c. As a hadron absorber it also 
detects showers with an angular resolution of a few degrees. 

The Belle detector has operated reliably over the first two years of KEK-B 
operation, accumulating 78 fb^^ on the T(4S') resonance corresponding to 85 
million BJ3 pairs that can be used to study CP violation in decay.! 




Figure 4: A schematic view the ingredients necessary to observe CP violation in 
B meson decay at an asymmetric e^e~ colhder. 

3 The Measurement of sin2(/)i 

There are three main ingredients that go into the measurement of a CP phase 
- for example sin20i in J/tpKg decay - at Belle. These are illustrated in fig. 
First, one must identify the decay vertices of the two B mesons. Second, since 
the final state under study is a CP eigenstate, we must tag the flavour of one B 
meson when it decays. Since the T(4S') decay products evolve in a coherent J = 1 
state, the decay of one of the B mesons into a B^ (in the example shown in fig. ^) 
projects the other B meson into a known B^ state. Finally, one must identify 
a sample of candidates which are CP eigenstates - for example J/ipK^ decays. 
I will briefly describe each of these aspects of the measurement in the following 
sections. 

3.1 B Meson Decay Vertexing 

Identifying and separating the two B decay vertices on an event-by-event basis in 
Belle is a subtle process. Despite the boost from the asymmetric beam energies 
the B mesons only travel a few hundred microns {'jct ^ 200/im) on average before 
they decay. The silicon vertex detector in Belle measures the decay vertices with 
a precision of about 100 fim. A J/ipKg decay vertex is typically determined 
with a precision of 75 fim, slightly better than the flavour tag B decay, which 
has a precision of 140 /im, because there are generally more reconstructed tracks 
attached to the former vertex. 
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Figure 5: Left: The distribution of decay time differences for and mesons 
measured by Belle. The open circles are the data, the solid line is the result of 
a maximum likelihood fit that includes a component proportional to the lifetime, 
and two components (dashed and dotted) that parametrise the detector resolution. 
Right: The distribution of decay time differences for B^ and B~ mesons. 



The precision of Belle's vertexing is demonstrated by its B lifetime measure- 
ments (shown in fig. |^) where one sees a clear difference between the detector 
resolution (dashed line) and the longer lived B meson decays (solid line through 
the data points). These fits give precise measurements of the charged and neutral 
B meson lifetimes and their ratio: t^/tq = 1.09±0.03.§ The detector resolution is 
understood out to ten B lifetimes and over three orders of magnitude in B decay 
rate. This is one of the crucial ingredients to measuring CP violation. 

3.2 Flavour Tagging in B'^ Meson Decay 

A second key ingredient is being able to tag the flavour of the B meson accompany- 
ing the CP eigenstate in the BB decay. This is done without fully reconstructing 
the opposite B - which would result in a significant loss in efficiency. Instead we 
look for characteristics of the other B decay that tag its flavour. For example, high 
momentum leptons indicate a direct semi-leptonic B decay, where the charge of 
the lepton is determined by the flavour of B. Lower momentum leptons arise from 
cascade semi-leptonic decays where the B meson first decays to a charmed meson 
and the latter decays semi-leptonically. This results in the opposite correlation 
between the lepton charge and B meson flavour. Such information is combined in 



Information of charged tracks 



Slow pion 



Lambda 



Kaon 



Track-level 
look-up tables 



Lepton 



Select track 

with 
largest "r" 



q.r 



Calculate 
combined "q.r" 



(q.r)K/A 



Select track 

with 
largest "r" 



q.r 



Event-level look-up table 



Flavor information "q" and "r" 



Figure 6: A schematic view of the Belle B flavour tagging algorithm. Events 
are classified by the following characteristics: charged leptons, kaons, lambdas 
and slow pions. This information is sifted through a set of look-up tables (based 
on Monte Carlo simulation) to identify and quantify the most reliable flavour 
tagging information in each event. The result of this procedure is the assignment 
of a "charge", q, and the reliability, r, of the flavour assignment. 



a set of look-up tables shown in fig. ^ that classify potential flavour information 
according to the sign of the b quark charge, q, and the reliability of the tagging 
information, r. 

While the look-up table information: q, r and their correlations - when more 
than one piece of tagging information is available in a single event - is extracted 
from a Monte Carlo simulation, the efficiency of the tagging algorithm is calibrated 
on a data control sample. Using a sample of B ^ D*lv decays - where one knows 
the flavour of the B meson from the charge of the lepton in the final state - the B 
mixing parameter is measured from the time-dependence of the observed decays. 
We use our fiavour tagging algorithm, described above, to classify events into six 
different ranges of r, shown in fig. |^. One sees that the observed amplitude of the 
mixing oscillation is much less in the sample tagged with lowest reliability (top 
left) while the amplitude of the oscillation (at At = 0) almost reaches 1 for those 
events that our algorithm tags with the highest reliability (bottom right). 

Figure ^ shows the correct tag probability from the six plots in fig. versus 
the average r for each sample. The correct tag probability, 1 — 2w;, is extracted 
from wi - the probability that a tag gives the wrong charge sign. In cases were 




Figure 7: B mixing measurements in self-tagged —>■ D*^l~v samples. This 
sample calibrates our tagging efficiency. The amplitude of the mixing observed, 
hence the efficiency of our tagger, depends on the reliability of the tag, which 
increases from top-left to bottom-right in this series of plots. 




Figure 8: The correct tag probability (1 — 2wi) versus tag reliability, r. The 
wrong-tag fraction {wi) is measured from the observed mixing in self-tagged 
samples (see fig. 0). The measured tagging probabilities (1 — 2wi) are used to 
weight the flavour determination in fits to CP eigenstates. 



wi approaches 0.5 (a 50-50 guess at the charge of the b quark) the probability 
of correctly tagging the flavour of the B meson goes to 0. Conversely, when wi 
goes to 0, one approaches a 100% correct tag probability. Figure || shows a strong 
correlation between our tagger's reliability, r, and the correct tag probability as 
measured in our control sample. These measured tag probabilities are used to 
weight the CP eigenstate decays in the flt to extract their decay asymmetry. 

With this algorithm we extract some tagging information from 99.5% of all B 
decay candidates in Belle. While some of the events have tags with low reliability, 
we measure an overall tagging efficiency of 28.8 ± 0.6%, corresponding to almost 
1/3 of our B sample being perfectly tagged. This represents an improvement 
from an effective tagging efficiency of 27.7 ± 1.2% in previous measurementsi 
arising from an increase to our low-momentum track reconstruction efficiency and 
an improved silicon alignment - allowing us to associate more tracks with the 
tagging B decay vertex. The reduced uncertainty on our tagging efficiency is the 
result of tripling of the size of the control data sample. 



Mode 


CP {^k) Candidates Purity (%) 


J/V^ir^(7r+7r-) 
J/^/'ir^(7r07rO) 

XciKl 

J/^K*\KItx'') 


-1 1116 98 
-1 162 82 
-1 172 93 
-1 67 96 
-1 122 68 
1 (81%) 89 92 
1 1230 63 


Total 


2958 



Table 1: Summary of the CP candidates, after flavour tagging and decay vertex- 
ing, used to measure sin20i. 

3.3 CP Eigenstate Event Samples 

From 78 fb~^ of data taken at the T(4S') resonance Belle has identified a sample of 
almost 3000 ccK decay candidates. These are summarised in table |I] where 

one sees that a little over 1/3 of the candidates come from J/ipK^ decays where 
the is reconstructed in a tt+tt^ final state. This golden mode also has the 
highest purity. We also have a number of other modes that are used to measure 
sin 201 with other cc resonances and other final states. Finally, about 40% 
of our sample is in the form of J/tpK^ decays that have the opposite CP eigen- 
value. This sample provides an important cross-check for possible CP dependent 
systematic effects in our measurement. Figure | shows the mass distribution for 
our candidates with candidates in the final state. 

While our K^^ detector does not measure the energy of the hadronic shower of 
a neutral Kaon decay with high precision, it does measure the shower direction 
with an accuracy of a few degrees. Having fully reconstructed a. J/ip meson one 
can hypothesize that a ^ J/ ipK^ decay has occurred and infer the magnitude 
of the momentum using only the direction, by constraining the two body 
system {J/il)., K^) to have the B^ mass. Imposing this constraint reduces the 
number of handles one has to reject background but still provides one degree 
of freedom - chosen to be the momentum of the B^ candidate in T(4S') centre 
of momentum (shown in fig. |TU|). Successfully reconstructed J/ipK^ candidates 
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Figure 9: The beam constrained mass for tlie different 
used to measure indirect CP violation in Belle. The largest sample comes from 
the J /'ipK^{TT^TT~) decays (dot-dashed line). The other modes combine to make 
up about 20% of our sample (see table p. 



peak near ^ 0.3 GeV/c, while backgrounds generate a flatter distribution in 
^CTns_ pj-om the distribution in fig. |10| we extract 1330 candidate events (before 
flavour tagging or decay vertex reconstruction) and estimate their purity to be 
63%. Since many of the backgrounds under the J/ipK^ peak result from other 
B meson decays care must be taken when extracting a CP asymmetry from this 
sample to account for the CP asymmetry of the background. We do this with 
Monte Carlo studies whose ability to constrain the CP content of the background 
is reflected in a systematic uncertainty on sin 20i from this channel. 
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Figure 10: The centre of mass momentum distribution for Belle's J/ipKi^ can- 
didates used to measure indirect CP violation. The signal forms a peak at 
p^*^ = 0.33 GeV/c, while backgrounds - where one or more final state particles 
have not been included in the candidate reconstruction - have their momentum 
smeared out. 




-8 -6 -4 -2 2 4 6 8 



At (ps) 

Figure 11: The time dependent asymmetry for and 5° decays to ccKg 
i^cp = ^1)- The amphtude of this oscillation, extracted from an event-by-event 
maximum likelihood fit, gives sin20i = 0.716 ± 0.083. 



3.4 Results 

We perform a maximum likelihood fit to all our 5° ccK candidates to extract 
their time dependent asymmetry (as shown in eqn. 0). This fit includes the flavour 
tagging probability, uncertainties on the measured separation between the two B 
meson decays (converted to At using the B meson boost) as well as the purity 
of the signal for each type of decay. The J/ipKg{n^7r^) candidates receive the 
highest weight because they have much smaller backgrounds (see table |I]). 

We first perform this fit separately for candidates with CP = —1 (the bulk 
of our data) and CP = +1 (our J/ipK^ sample), shown in figs. |TT| and |T^, re- 
spectively. We find the magnitude of the fit asymmetries are equal, but they have 
opposite sign - as expected. The fit results are reported in table |[ We then 
combine all candidates into a single fit (see fig. |T3|) - inverting the sign of At for 
the CP = +1 candidates - and obtain: 

sin 201 = 0.719 ± 0.074(stat.) ± 0.035(s?/s.) 

The systematic uncertainties involved in this measurement are listed in table 0. 
The largest of these uncertainties are derived from measurements of our control 



samples - for example the flavour tagging discussed in section 3.2 - and thus can 



be expected to shrink as more data becomes available.! The overall systematic 
uncertainty is small compared to the statistical precision which bodes well for 
future measurements. 
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Figure 12: The time dependent asymmetry for and decays to J/ipK^ 
(^cp = !)• The amphtude of this oscillation, extracted from an event-by-event 
maximum likelihood fit, gives sin20x — 

0.781 ±0.167. 
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Figure 13: The combined time dependent asymmetry for B^ and B^ decays to all 
ccK^ final states (flipping the sign of the candidates with ^cp = !)• An event-by- 
event maximum likelihood fit gives sin20i = 0.719 ± 0.074. 



Uncertainty Source Value 



Vertexing Reconstruction 0.022 

Flavour Tagging 0.015 

Vertex Resolution 0.014 

Fit parametrisation 0.011 

J/ipKl Background 0.010 

Amd, tb < 0.010 

Total 0.035 



Table 2: Systematic uncertainties on the measurement of sin20i. 



3.5 Cross-checks of sin20i 

We perform a number of cross-checks on non-CP eigenstate samples, all of which 
exhibit null asymmetries with statistical precisions ranging from 0.02 to 0.09. We 
do not include these as systematic uncertainties as we find no evidence of bias 
but they provide additional confidence in our main result. We have also sub- 
divided our data into different sub-samples to see whether there is any evidence 
for systematic variations in our result. Finding none (see table ^ gives us fur- 
ther confidence in our main result but again do not ascribe additional systematic 
uncertainty as the variations are all consistent with our quoted result within the 
statistics of the smaller sub-samples involved. 



4 CP Asymmetries in B ti^ti Decay 

Having observed indirect CP violation in B ^ J/ipK decays it is natural to 
study other CP eigenstates. The decay tt^-k^ is interesting because there 

are at least two significant amplitudes that can interfere in the direct decay. These 
amplitudes are show diagrammatically in fig. |T^. With more than one amplitude in 
the direct decay path one can have a cos AmAt dependence in the CP asymmetry. 
This dependence can come from the interference of the direct amplitudes and can 
appear in addition to the sinAmAt time-dependence (compare for example to 
eqn. |1|) from the interference with the mixing amplitude: 



Subsample (stat error only) 



J/ipKli'n+'n-) 0.73 ±0.10 

{cc)K^s (except J /i)Kl{'K+^~)) 0.67 ± 0.17 

J/i^Kl 0.78 ±0.17 

ftag = B'' 0.65 ±0.12 

ftag = 0-77 ±0.09 

r < 0.5 1.26 ±0.36 

0.5 <r< 0.75 0.62 ±0.15 

0.75 <r 0.72 ±0.09 

All 0.72 ±0.07 



Table 3: Cross-checks of sin 20i from fits to subsets of the data. 



^(B°^Fcp)-f (B°^Fcp) 
^(B«-Fcp)+^(B0~>^'CP)' 



Acp{At) = 

= Sf sin AmAt + Apcos AniAt. [2] 



The K^TT~ branching fraction - three times larger than the 7i^tt~ 

branching fraction - is strong evidence that the penguin amplitude (right diagram 
in fig. |l^) is not small. In the decay B^ tt+tt^ the CKM elements involved 
predict an asymmetry proportional to the angle 02- However the penguin diagram 
can introduce a non-CKM phase thus: 

Sf ~ sin (202 ± 0) 

where 6 can come from phases at the gluon vertices in the penguin process. While 
waiting for measurements of the various B ^ nir branching fractions! with suffi- 
cient precision to constrain 6, it is still interesting to measure the tt+tt^ asymmetry 
and perhaps get some hint for the size of the combination sin(202 ± 0)^ 



4.1 TT+TT Candidate Selection 

The candidates and result presented here are from 43 fb~^ of data collected by 
Belle.! While the B vertexing and fiavour tagging described in sections 
can be used in a measurement of the vr+vr" asymmetry, the candidate selection 
is more challenging. The branching fractions for B h^h~ are an order of 
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Figure 14: Two tree level diagrams for tt+tt" decays. Given that there are 

two paths to the same CP eigenstate, their interference can produce direct CP 
violation in this mode. 



magnitude smaller than those of the ccK final states - on the order of 10"^. 
Furthermore the B — > h'^h^ decays have only two, relatively high momentum, 
tracks in the final state. Such candidates are much more prone to look like the 
e'^e" qq continuum background. Belle uses a set of event-shape variables 
(Fox- Wolfram moments) and the reconstructed B^ candidate direction to suppress 
continuum background. These are combined into a single variable, R, shown in 



figure |T5| for BB data (open circles) and off-resonance data (closed circles). By 
requiring R > 0.8 we reduce the continuum background by an order of magnitude 
while retaining two thirds of our B candidates. 

Since these events are fully reconstructed we are able to measure both the 
beam-constrained mass, mtc, and the difference between the reconstructed energy 
of the B candidate and the beam energy, AE. Figure |TB| shows the distributions 
of both of these variables for our B h^h^ sample. Comparing the mbc distri- 
bution for these candidates to that of the ccK candidates in fig. ^ one sees that 
the background is much more prominent. In the AE distribution (on the right in 



fig. |T6D one can see the source of these backgrounds. The remaining continuum 
background populates AE almost uniformly - falling linearly with increasing AE 
due to phase space - and is represented by the dotted line. B decays to three 
or more particles, where we have missed one in our reconstruction, populate the 
low AE region (grey line) at a much lower rate than the continuum. Finally, 
B^ K^-K^ candidates that have not been rejected by our particle identification 
produce a peak at AE ^ —0.045 GeV. This peak is shifted because these candi- 
dates have been reconstructed assuming that both B^ decay daughter tracks are 
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Figure 15: Low multiplicity final states are subject to higher background from 
the e~^e~ continuum. To suppress these Belle combines event shape variables to 
get an overall likelihood, R. The solid points are continuum data, while the open 
points are a control sample of fully reconstructed, low-multiplicity, B candidates. 
— > 7r+7r~ candidates are required to have R > 0.8. 

pions (with a mass of 139 MeV/c^). If one of them is really a kaon (with a mass of 
494 McV/c^) then the reconstructed energy shifts down by about 45 MeV. This 
background is shown by the dot-dashed line. 

Wc restrict our CP fit to candidates that have \ AE\ < 0.067 GcV, eliminating 
the mis-rcconstructcd B background entirely and leaving us with 74 7r"'"7r~ candi- 
dates, 28 K^Ti^ candidates and about 100 candidates that come from the e^e" 
continuum. 

4.2 The CP Fit 

We perform an event-by-event likelihood fit to the time dependent asymmetry 
for the B — > tt'^tt" candidates. In order to maximise our sensitivity to the tt'^tt" 
asymmetry this fit weights the events according to their continuum likelihood 
(events near R — 1 are weighed more heavily) as well as their fiavour tagging 
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Figure 16: Left: The beam constrained mass, mtc, for the tt+tt^ candidates 

used to extract the CP asymmetry. Right: The energy difference distribution, 
AE, for the same candidates. In AE the 7r"*'7r~ signal peaks at (dashed Gaus- 
sian), the K~^'n'~ background peaks at -0.045 GeV (dot-dashed Gaussian). The 
continuum is parametrised as a falhng hnear background (dotted hne) while 
decays with more that two particles in the final state appear at low AE. 



probability and vertex reconstruction quality - as was done for sin20i. 

Figures |1^ a) and b) show the At distributions for events tagged as coming 
from B^ and B^ decays, respectively. The continuum and K^tt' backgrounds 
are represented by the dashed lines, symmetric about At = in each case. The 
hatched areas in each plot are what the fit ascribes to the B tt+tt^ candidates. 



The sum of these two distributions give the solid line which fits the data. Figure |T7 
c) shows the At distribution for just the B tt+tt^ component of the fit. Here we 
see that our fit predicts that there are almost twice as many B^ as B^ candidates 
in our sample - independent of At - an indication of direct CP violation. Finally, 
fig. d) shows the asymmetry between the two distributions in fig. ^ c). Here 
one also sees a At dependent asymmetry leading to our fit result: 

=-1.2ll°:f7(stat)l°:l^(sys), 
A^^= 0.94j:[Ji^(stat)±0.09(sys). 

This is three sigma evidence for both direct CP violation in the decay B ti^ti~ 
and indirect CP violation.! 

We have performed a large number of cross-checks on this fit, including fits 



to the K^TX^ sample which should not have a CP violating asymmetry 

and fits to continuum side-bands (fig. |1^ e). None of these fits show significant 
asymmetries. We have simulated the size of our statistical uncertainties. We 
find that our statistical uncertainty on S^^.,^ is somewhat smaller than might be 
expected for a sample this size. Still there is a 5% chance that we could have 
gotten a smaller uncertainty. While unlikely, this is not an unacceptable statistical 
fluctuation. Finally, we have performed a series of toy Monte Carlo simulations, 
each generated with (5,^77, A-k-k) = (0, 0) and sample sizes (signal and background) 
that mimic those in our data. We find that 1.6% of these toy Monte Carlo samples 
return a fit value farther from {St^t^, Aj^j^) = (0, 0) than our data. It will be 
interesting to see how this result evolves as more data becomes available. 



5 CP Asymmetries in Rarer B Decay Modes 

Having looked at two of the more plentiful B meson decays that are CP eigen- 
states, we turn to the next most likely B decay modes to show CP phases. One 
such class of decays is shown in fig. 0. With no tree level contribution these 
decay processes might be expected to have branching fractions smaller than those 
discussed above. In fact, the mode B —>■ rj'K^ has a branching ratio that is some- 
what larger than the mode B tc^tc" and clearly larger than the, as yet unseen, 
decay B r]K. That the rj'K branching ratio is so large (about 5 x 10~^) is 
not easy to understand theoretically - hinting that there may be more that just 
a charged W boson involved in the loop at the top of the diagram. If there were 
contributions from a charged Higgs to this amplitude, it could introduce an addi- 
tional phase, interfering with the weak phase predicted by the CKM model. In the 
absence of such additional phases this decay would provide another way - albeit 
a much lower statistics way - to measure sin20i. Thus a comparison of the CP 
asymmetry in the decay f]'Kg to that seen in ccKg could help unravel the mystery 
of the large branching of B mesons into this mode and may provide a glimpse of 
physics beyond the Standard Model. 

5.1 CP Asymmetry in the Decay I'Kg 

Figure |TU] shows the beam constrained mass distribution for 128 candidate B^ 
rj'Kg decays from 78 fb^^ of Belle data. Here the t]' candidates have been re- 



Mode 


(sin20i)cfT Candidates 




0.76 ± 0.36(stat) ± 0.06(sys) 128 
-0.73 ± 0.64(stat) ± 0.18(sys) 35 
0.52 ± 0.46(stat) ± 0.11(sys)t[}:g^(CP) 95 



Table 4: CP fit parameters for b — »• sss decay modes, 
constructed in rj tt+tt^ and intermediate states and then combined with 



—>■ TT+TT combinations. The plot on the right of fig. [T^ shows the B^B^ 



asymmetry, as a function of At, where the candidates have been flavour tagged 
and vertexed in the same way as our ccK^ measurement of sin20i. The coeffi- 
cient for indirect CP violation derived from this fit (eqn. ^ with A = 1) is listed 
in table |[ While the statistics are still poor, this is a demonstration that we can 
measure CP asymmetries in these rare modes. More details on this analysis can 
be found in reference!! that describes this measurement made on the first half of 
the data shown here. 

We perform the same analysis on samples of B^ K^K^Kg decays. Though 
the number of candidates in this sample is even smaller, we can begin to probe 
their CP asymmetry. We separate the K^K~ candidates into two samples. When 
fTT-K+K- ~ fTT'cp the final state has CP = — 1 providing another potentially clean 
determination of sin20i. The statistics limited result is shown in table ^. The 
remaining sample is somewhat larger, providing a more incisive measurement of 
the CP asymmetry. Studies of similar non-resonant decays related to this one by 
isospin (such as B^ K^K^K^) have shown that our (ir+ir~)non-res-K^5 sample 
is 97lfg% CP even. This bound on the CP symmetry of the non-resonant K^K~ 
final state translates into an additional systematic on the effective sin 20, listed in 
table ^ The measurement of the CP eigenvalue in this mode will improve with 
additional data so the non-resonant mode might eventually provide another clean 
measurement of sin20i, or a glimpse of a non-Standard Model phase this decay. 



6 Summary and Future Prospects 

Belle has studied the decay of B^ mesons into a number oi CP eigenstates. The 
most precise measurement of sin20i to date comes from the modes B^ ccK. 



Most of the systematic uncertainties involved in that measurement come from 
cahbrations of our flavour tagging and decay vertexing procedure based on control 
samples. Thus it is reasonable to expect that the precision on sin20i - both 
statistical and systematic - will improve as more data becomes available. Figure pO| 
shows a prediction of how the precision on sin 20i will scale as Belle collects more 
data. The lowest curve shows the evolution of the systematic uncertainty on 
sin 201 in the ccK measurement. It should improve until as least 1000 fb~^ (1 
ab~^) of data has been collected. At that point we expect an overall precision 
(uppermost thin line in fig. 

The rarer modes will also benefit from the increase in statistics. Belle expects 
to be able to accumulate 1 ab~^ of data in the next 3 or 4 years. At that point the 
precision on the Tj'K^ asymmetry should approach the precision that we currently 
have on the ccK asymmetry. If there are other modes contributing to the phase 
in Tj'K^ decay un-natural fine-tuning would be necessary for our precision not to 
reveal some discrepancy between what would otherwise be two measurements of 
the same quantity. 

The B factories were built to confirm the CP violation predicted by the CKM 
model in meson decay - they have now achieved that milestone. Having made 
a measurement of sin 201 with a precision of 10% in the ccK^g modes, and are 
now pushing to improve the precision to a few percent. At the same time they 
are expanding the scope of their study to other decay processes. In doing so they 
will measure other angles of the unitarity triangle. This will over-constrain the 
CKM model and test whether there is physics beyond the three known families 
of quarks and their Standard Model electroweak decays. In the coming years the 
flavour sector will become ever more constrained reducing the range of physics 
that could explain why we have three families of quarks. 
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Figure 17: The results of the vr+vr^ CP asyiiimetry fit. a) The At distribution for 
candidates tagged with q = 1. The signal is shown as the hatched region, while 
the background (dominated by the continuum) is shown by the dashed line. The 
solid line is a projection of the overall fit and the points are the data, b) The At 
distribution for candidates with tagged q = —1. c) The circles (triangles) are the 
background subtracted data for q = (— )1 while the curves are fit projections for 
the two flavours, d) The difference between the two data-samples in c). e) The 
asymmetry observed in a continuum sample. 




q = u, c, t 



Figure 18: The Feynman-diagram for the decay v'^s- The unexpectedly 

large branching fraction into this mode hints that there may be extra contributions 
to the amplitude - which could introduce additional CP phases. 




Figure 19: Left: The beam constrained mass distribution for —>■ rj'Kg candi- 
dates. Right: The CP asymmetry as a function of At for these candidates. 




Figure 20: Prospects for improving the precision on sin 201 from the decay ccK^ 
(lower set of thin curves) with additional statistics. For comparison the precision 
expected on rarer modes {ti'K^ thick grey curve and (pK^ thick black curve) are 
also shown. 



